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Abstract. We study sample path deviations of the Wiener process from three 
different representations of its bridge: anticipative version, integral representation 
and space-time transform. Although these representations of the Wiener bridge 
are equal in law, their sample path behavior is quite different. Our results nicely 
demonstrate this fact. We calculate and compare the expected absolute, quadratic 
and conditional quadratic path deviations of the different representations of the 
Wiener bridge from the original Wiener process. It is further shown that the pre- 
sented qualitative behavior of sample path deviations is not restricted only to the 
Wiener process and its bridges. Sample path deviations of the Ornstein-Uhlenbeck 
process from its bridge versions are also considered and we give some quantitative 
answers also in this case. 



Let (Wt)t>o be a standard one-dimensional Wiener process on a filtered probability 
space (fl, J 7 , (J r t)t>o, -P); where the filtration {Tt)t>o is the usual augmentation of the 
natural filtration of the Wiener process W (see, e.g., Karatzas and Shreve [T8l Section 
5.2.A]). We consider the following versions of the Wiener bridge from a to b over the 
time-interval [0,T], where a,b G R (see, e.g., Karatzas and Shreve [T%j Section 5.6.B]): 
1. Anticipative version 
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1. Introduction 




< t < T. 



l 
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2. Integral representation 



a+(b-a)^ + £L-±dW 8 ifO<t<T, 



b if t = T. 

3. Space-time transform 

a + (b-a)l= + ^^Wjt_ if 0<t<T, 



if t = T 



Here the attribute anticipative indicates that for the definition of W* v we use the 
random variable Wt, where the time point T follows the time point t. In the 
sequel we will use the notation (W t br ) t6 [ 0) T] if the version of the bridge is not specified. 
All the bridge versions above are Gauss processes with the same finite-dimensional 
distributions. This can be easily calculated, since the versions all have mean function 
E(W t br ) = a + (b - a)±, < t < T and for < s < t < T we have the covariance 
function 

Cov(W: v , W t av ) = Cov (w s - ^W T , W t - !=w T 



st st st T — t 
^ T T T ^ T * 



C0V(W:, Wf) = COV ^-i dW r , £ ^ dW^j 



(T-s)(T-t), /rjl ,/ 1 1\ T-t 

Ur = (T-s)(T-t)[—---)=s — 



(T-r) 2 v /v \T — s T 



and 



(T-s)(T-t) sT T- 

— s 



T 2 T-s T ' 

where we used that the function [0,T) 3 t i— >■ is monotone increasing. Altogether, 
for all < s < t < T we have 

(1.1) Cov(W?,W?) = s^. 

We note that the finite dimensional distributions of the above Wiener bridge ver- 
sions coincide with the conditional finite dimensional distributions of the Wiener 
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process (a + W^)te[o,T] starting in a and conditioned on {Wt = b}; see, e.g., Problem 
5.6.13 in Karatzas and Shreve |18j or Chapter IV.4 in Borodin and Salminen [7J. 
Bridges of Gaussian processes have been generally defined by Gasbarra et al. |14| . 
while from the Markovian point of view the reader may consult Fitzsimmons et al. 
|13| . Barczy and Pap [I], Chaumont and Uribe Bravo [9j, and the more recent Bryc 
and Wesolowski [8J which deals with the inhomogeneous case. 

Moreover, it follows from the definitions that all bridge versions have almost sure 
continuous sample paths. The (left) continuity of the trajectories at t = T is not 
obvious in case of the integral representation and space-time transform. Corollary 
5.6.10 in Karatzas and Shreve |18) yields the desired continuity for the integral repre- 
sentation, whereas the strong law of large numbers for a standard Wiener process (see, 
e.g., Problem 2.9.3 in Karatzas and Shreve [IB]) for the space-time transform. Hence 
the anticipative version W av , the integral representation W a and the space-time 
transform W st induce the same probability measure on (C[0, T], B(C[0, T])), where 
C[0,T] is the space of continuous functions from [0, T] into K and B(C[0,T]) 
denotes the Borel a-algebra on C[0, T\. This underlines and explains the definition 
of a Wiener bridge from a to b over the time-interval [0, T] (see, e.g., Karatzas 
and Shreve [LSI Definition 5.6.12]), namely, it is any almost surely continuous Gauss 
process having mean function a+ (6 — a)^, t G [0, T], and covariance function given 
in flTTj ). 

Furthermore, according to Section 5.6.B in Karatzas and Shreve [18j or Example 
8.5 in Chapter IV in Ikeda and Watanabe [16], the above versions of the Wiener 
bridge are solutions to the linear stochastic differential equation (SDE) 

h — U/br 

(1.2) dWf T = — — l — dt + dW t , < t < T , with W$ r = a. 

By Theorem 5.2.1 in Cksendal |19j or Theorem 2.32 in Chapter III in Jacod and 



Shiryaev [22], strong uniqueness holds for the SDE (1.2), and (W")^^^) is the 
unique strong solution of this SDE being adapted to the filtration (J r t)t&[o,T)- Whereas 
(W / i av ) tg [ ,T) is only a weak solution to the SDE (1.2); it can not be a strong so- 
lution, since the definition of the anticipative representation formally requires in- 
formation about Wt, although Wf and Wt are independent for every t G [0,T] 
(indeed, Cov(W / t av , Wt) = Cav(W t , W T ) - | Cov(W T} W T ) = 0, t e [0,T]). The 
space-time transform representation {Wf) t& [o^T) is only a weak solution to the SDE 
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(1.2), too, since it is adapted only to the nitration (J 7 tr )te[o,T) and T tr ^ JFt, 
t G (0,T). We also note that, even though the three bridge versions have the same 
law on (C[0, T], B(C[0, T])), their joint laws together with the Wiener process through 



which they are constructed, are different (see Propositions 2.1 and 2.4). Our aim is 
to elucidate their sample path deviations compared to the original Wiener process 
(a + Wt)te[o,T\ starting in a. A motivation for our study is given at the end of this 
section. 

According to simulation studies, for a typical sample path of the Wiener process the 
deviations from its anticipative bridge version and its space-time transform are larger 
than from its integral representation of the bridge; see Figure [Tj Note that in general 




FIGURE 1. Two typical sample paths of the Wiener process (rows, thick 
lines) and its deviations from the anticipative version (left column), the in- 
tegral representation (middle column), and the space-time transform (right 
column) of the Wiener bridge from to over the time-interval [0, 1]. 



the deviation from the space-time transform bridge version is hard to compare with 
the other deviations, since (Wf)te[T/2,T) depends on the non-visible part (Wt)t&[T,oo) 
of the Wiener process. Our aim is to give quantitative answers to this qualitative 
behavior observed from simulation studies and thus to study the path deviations on 
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[0,T): 



1.3) 



a + W t - Wf 



a + W t - Wf 



(a-b)^ + ^W T , 



a + Wt-W, 



(« - 6) j; + 
- b) - + 



t 



T-s 
T-t 



T 



WjT_ 
T-t 



Note that the dependence of the path deviations in (1.3) upon the starting and 



endpoint of the bridge (a and b) is only via their difference a — b. Hence without loss 
of generality we can and will assume a = in the sequel. 

Simulation studies also show that the above typical behavior is reversed in case 
the endpoint Wt of the Wiener sample path is close to the prescribed endpoint b of 
its bridge, namely, for such a sample path of the Wiener process the deviation from 
its anticipative bridge version is smaller than from its integral representation of the 
bridge or from its space-time bridge version; see Figure [2j We aim to give quantitative 




FIGURE 2. A sample path of the Wiener process with W\ ~ (thick line) 
and its deviations from the anticipative version (left), the integral represen- 
tation (middle), and the space-time transform (right) of the Wiener bridge 
from to over the time-interval [0, 1]. 



answers to this effect and thus in Section 2 we will particularly compare the so-called 
expected p-th order sample path deviations 



3j J \W t -W^\ p dt\ = J E(\W t -W^\ p )dt 
for p = 1, 2 and in case of p = 2 we will explicitly calculate the conditional analogue 



E 



J {W t - Wf 1 ) 2 dt W T = d^j = E({W t - W^ T f | W T = d) dt 
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for prescribed endpoints Wt = d, d £ R of the original Wiener process. In the 
above formulas, integration over the time-interval [0, T] and taking expectations can 
be interchanged. Indeed, since we have continuous sample paths, we can consider 
monotone approximations of the integrals by Riemannian sums with nonnegative 
summands and then apply the monotone convergence theorem for conditional expec- 
tations. In what follows expected first and second order sample path deviations will 
be called expected absolute and quadratic path deviations, respectively. 

We will further show in Section 3 that the above mentioned qualitative behav- 
ior of sample path deviations is not restricted only to the Wiener process and its 
bridge versions: sample path deviations of the Ornstein-Uhlenbeck process from its 
bridge versions are also considered. Here we give some quantitative answers, too, see 
Theorem 13.71 

In the Appendix we present an auxiliary result which is used for proving almost 
sure continuity of the integral representation of the Ornstein-Uhlenbeck bridge at the 
endpoint of the bridge. 

Our results are to be seen as paradigmatic examples that give rise for future work 
concerning more broad questions of how certain pathwise constructions of Gaussian 
or Markovian bridges can differ, although they obey the same law. The reason for 
concentrating on the Wiener and on the Ornstein-Uhlenbeck process here is the pos- 
sibility of giving explicit expressions for some quantities (such as second moment) re- 
lated to the path deviations of different bridge versions to the original process through 
which they are constructed. In particular, the case of an Ornstein-Uhlenbeck process 
shows that explicit expressions for path deviations can soon become unwieldily. As 
a future task, one may also address the question of existence of a bridge version that 
minimizes the distance to the unconditioned stochastic process in a certain sense. 
Moreover, one may present other indicators for different sample path behavior of 
the Wiener and Ornstein-Uhlenbeck bridge versions, such as Hellinger distance, and 
address the question for more general process bridges. 

To further motivate our study, we point out that similar problems were considered 
by DasGupta |12| . Bharath and Dey [2j and Balabdaoui and Pitman 1 1 . Namely, 
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DasGupta (TJJ Theorem 1] gave an infinite series representation of the expectations 

E Qf |Wj br - fd - W t \ dt^j , 5 6(0,1], fieR, 

where (Wt)t e [o ii and (W} >r )te[o,i] denote respectively a standard Wiener process and 
an independent Wiener bridge with a = b = and T — 1. For some special values 
of 5 and // the exact values were also calculated. The motivation of DasGupta for 
calculating the expectations above is to understand whether distinguishing between a 
Wiener bridge and an independent Wiener process with possible drift on the basis of 
observations at discrete times is intrinsically difficult. It turned out that distinguish- 
ing one from the other is not an easy task. DasGupta studied the likelihood ratio test 
for testing the null-hypothesis H : X t = W^, t G [0, 1], against the alternative hy- 
pothesis H\ : X t = Wt+nt, t G [0, 1] for some /j G R, based on discrete observations 
from a process (X t )t^[o,i}- Recently, the question of distinguishing a Wiener process 
from a Wiener bridge was also considered by Bharath and Dey [2J. Note that in our 
setup (Wt)t 6 [o,i] and (W f br )t e [o,i] are not independent. Hence our results may be useful 
to answer the question of distinction in case the Wiener bridge is constructed by the 
help of the original Wiener process and not an independent copy. One can address 
the same question for Ornstein-Uhlenbeck bridges or for more general process bridges. 
Our calculations in the Ornstein-Uhlenbeck case can be considered as a first step to- 
wards the corresponding calculations of Section 2 in DasGupta |12| . Balabdaoui and 
Pitman |T] gave a representation of the maximal difference between a Wiener bridge 
and its (least) concave majorant on the unit interval. As an application, expressions 
for the distribution, density function and moments of this difference were derived. 

The presented results might also be applied to the study of animal movements. 
Home et al. |15j use a two-dimensional Wiener bridge to model the unknown 
movement of an animal between two consecutively observed positions of the ani- 
mal. The model is used to investigate questions on the mean occupation frequency 
E(i J o r l A {X\% X^ t ) dt) in a region A G B(R 2 ), where (X^) te[0jT] and (X% t ) te[0 , T] are 
independent Wiener bridges such that (Xj^X^o) and {X\ X T , X^j) are the starting 
and ending positions of the animal at time and T, respectively. If the region A 
depends on the original (independent) Wiener processes (Xi^)te[o,T], (X2,t)te{o,T\, e.g., 
for questions concerning the closeness of the animal's path to the path of a Wiener 
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process, our results show that the expected occupation frequency heavily depends on 
the chosen version of the bridge. 



2. Path deviation of the Wiener process from its bridges 

2.1. An indicator for different sample path behavior of Wiener bridge ver- 
sions. A first indicator for different sample path behavior of the bridge versions is 
the correlation function g(W^ T , W t ) of these bridge versions and the original Wiener 
process. Note that {W^, Wt)te[o,T] is a two-dimensional Gauss process and the 
correlation coefficient of the two coordinates is given in the next proposition. 

Proposition 2.1. For all t E (0,T), we have 



fj{] r;>\\\y j = o(W? : W t ) = \[^ and g (wi\w t ) = \og ' 



T t °T-t 



Proof. By (1.1), we get for every < t < T 

Var(^ br ) = Cov(^ br , W?) = t^. 
We easily calculate for every < t < T 

Cov(WT, W t ) = Cov(W t> W t ) - % - Cov(iy T , W t )=t-^= t 7 ^, 
Cov(W7, W t ) = Cov (J* ^ dW s , J* 1 dW s ) = ^ ^ ds = (T — t) log ^- 
and 

Cov(Wf , W t ) = Cov H^-W^, WA = t 7 ^. 
Thus we get for every < t < T, 



1 2/1 ) g(Wr, W t ) = , f ~ = J 7 ^- = 6 {W? , W t ) 

(T-t)-t V 1 



and 



% (T-t)-t 1 T ~ t 

concluding the proof. □ 
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Remark 2.2. For all T G (0,oo), the function (0,T) 3 t \-> g(W^,W t ) is strictly 
decreasing. For the anticipative version and space-time transform, it is an immediate 



consequence of (2.1). For the integral representation, it is enough to check that 

dt { 1 l0g T—t)= * < ° 

for all t G (0, T), which is equivalent to show that 

hit) := [T - log + * < 0, t G (0, T). 

Using that log(l — x) = — J^jfcLi t~ ^ or a ^ — 1 < x < 1, we get 

k oo k 



v / fc=i v / fc=i fc=i 

^ V(A; + l)T fc 2kT k )~ \k + 1 2A? J 



fc=i v v ' 7 fc=i 

Note also that £(W br , W t ) -»• 1 as t | 0, and 0(W t br , W t ) -» as t t T. Hence W f bl 
and Wt, £ G (0,T), are positively correlated for all bridge versions. Moreover, 



/ \ JT(T - t) T lT-t 

(2-2) ^ >-\og—>J — , tG(0,T). 



Indeed, (2.2) is equivalent to — |- > log (l — |;J for all t G (0, T), which follows by 
log(l — x) < —x for all < x < 1. 

Hence the integral representation is more positively correlated to the original process 
than the anticipative version and the space-time transform. □ 

2.2. Gauss and conditional Gauss distribution of path deviations. First we 
study the distribution of the path deviation W t — W t , t G [0,T). 



Proposition 2.3. Let (W t )te[o,T] be a Wiener bridge from to b over the time- 
interval [0,T], where 6 el. Then for all t G [0, T), the path deviation W t — W^ r 
is normally distributed with mean E(W t — W^ r ) = —bk and with variance 



Vax(Wi - WD = Var(H/ i - W?) = %, 
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Var(^ t - Wn =t(l + + 2(T - t) log ^* =: a 2 (t). 

Proof. With a = 0, by (1.3), for every < t < T the path deviation Wt — W^ r is 
normally distributed with mean E,(W t — W^ r ) — — and with variance 

T" x ) T 7 
* ' T-t\ 2 , /•* T-t fT-t\ 2 



Var(W( - W t av ) = Var ( ^W T ^ 



t + 2(T - t) log ^ + (T - i) 2 - I ) = a 2 (t), 



and 



Var(W t - Wf ) = Var - ^^W^J = Var - W t ) + 

T-t\ 2 (tT \ t 3 _ (T -t)t 2 + t 3 _t 2 



T ) \T — t I T 2 T 2 T 

concluding the proof. □ 



By Proposition 2J3, for every < t < T, the variance of the path deviation of the 
integral representation from the original Wiener process is smaller than those of the 
anticipative version or the space-time transform, since we have a 2 (t) = 2t-^ + 2(T- 
t) log(l — |;) and thus 

(2.3) o- 2 (t)<^ te(0,T). 



Indeed, (2.3) is equivalent to — ^ > log (l — 4) for all t G (0, T), which holds, since 
log(l — x) < —x for all < x < 1. 

Next we examine the conditional distribution of the path deviation Wt — W^ r 
given the endpoint Wt- 

Proposition 2.4. Let (W^ T )t£\o,T] be a Wiener bridge from to b over the time- 
interval [0, T], where 6g1. Then for all te[0, T) and d G M., the conditional 
distribution of the path deviation Wt — W^ r given Wt = d is normal with mean 

(2.4) E(W t - W™ \W T = d) = {d- b)^=, 
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(2.5) E(W t - W* | W T = d) = {d-b) t - + log 

t d 



(2.6) E(W t - Wf | W T = d) = -b- + -(2t - T) ■ l w (t), 
and un£/i variance 

(2.7) Var(W t - W t av | W T = d) = 0, 

T — t (T -t) 2 T -t 

(2.8) Var(W^ - | W T = d) = 2t— + log ^ 

(T-tf f T-O 2 
log 



T V T 

f2 (Of _ Tl 2 

(2.9) Var(Wi - W? \ W T = d) = - - 1 - ) ■ l |?fl (t). 

Proof. For all < t < T, the joint distribution (W t — W^ 1 , Wt) of the path devia- 
tion and the endpoint is a two-dimensional normal distribution and, by Theorem 2 
and Problem 5 in Chapter II, §13 of Shiryaev |22j . it is known that the conditional 
distribution of Wt — W^ r given Wt = d is normal with mean 

(2.10) E (W t - W?) + d ~®$ff Cov (W t - < r , W T ) 
and with variance 

^ (Cov (W t - Wf\ W T )) 2 

(2.H) v *(w,-wn- { 4^ T " ■ 

Here we have 

Cov(iy 4 - W? v , W T ) = Cov (^W T , W T J = t, t G [0, T), 



and thus (2.10), (2.11) and Proposition 2.3 yield that 

E(W t - Wr I W T = d) = -b l - + ±t=(d- b)^ 

Vav(W t - Wt v \W T = d) = ^-^=0. 
We note that the above formulae follow immediately, since in case of Wt = d, we 



have W t — W^ v = (d — b)^. Further, we have 



Cov(W t - W'r W T ) = Cov QT L^- dW s , £ 1 dW s \ = £ ^-i ds 

* / T — t\ T — t 

1-^ d S = t + (T-t)log— - , 



o 



T — s 



12 
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and thus (2.10), (2.11) and Proposition 2.3 yield that 

d 



E(W t - W* | W T = d) = -b l - + r 



t+(T-t) log 



T-t 
T 



T-t\ 2 



Var{W t - | W T — d) — a 2 (t) 



(t + (T-t)bg^) 
T 

t 2 T-t 
2t-- + 2(T-t)log — 

---2(T-t)-log— - 



(T-f) s 



log 



T-t 
T 



This implies ( |2.5 ) and (2.8). Finally, we have 

Cov(Wi - Wf, W T ) = Cov(Wt, W T ) 

T-t 



mm 



T 

T-t tT 
T T-t 

T-t rp _ 

T 



T — t 

T 
tT 

T-t 1 




Cov (Wrr , W T 

\ T-t 



T 



if < t < f , 
2t — T if | < t < T, 



and thus (2.10), (2.11) and Proposition 2.3 yield (2.6) and (2.9). 



□ 



2.3. Comparison of tail functions. The tail of a normally distributed random 
variable Y„ a 2 with mean and with variance a 2 > has the form 



T llja2 (x)=F{\Y^\>x} 
= 1 -P 



1 



x — fi Y^ a 2 - fi x - (j. 



(7 

/! + X 



a 



(7 

fl — X 

a 



a 



1 - $ 



x — fx 



a 



x + /i 



(7 



x > 0, 



where $ denotes the standard normal distribution function. Since, by Proposition 
2~3~[ E(Wi - iy 4 br ) = t G [0,T), if we want to use the monotonicity in 

to show different behavior of the tails of the deviations Wt — W^ r , then this tail 



function should be an increasing function in a > for every fixed x > and fixed 
/i := -b± e R. We have 



T 
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2 1 / la; 2 + u 2 \ fx . fxa\ u . fxa 

In case < x < fi ^ 0, we have 

X /Xfi\ fi (Xjl\ X — H fXUL\ x + jl ( x\i\ 

a COSh UJ - a I W = ~2T 6XP I W + 6XP b) 

as a | 0. This shows that in general the tail function T„ 3.2(2;) is not increasing in 

cr > and thus is in general not helpful to analyze the different behavior of path 

deviations. 

In special situations such as b = = /i it is evident that (0, oo) 3 a i— > 7 0a 2^x) is 
strictly increasing for every x > 0. In this special case it follows immediately from 



the formula E(|y 0)(T 2| p ) = f^°x p 1 T (J 2(x) dx, p > 1, and (2.3) that for every p > 1 
and < t < T we have 

E (\W t - Wf\ p ) < E (I Wt - W t av | p ) = E (|Wt - . 
As a further consequence we get for every p > 1 

e ^ jT |w t - W7i p d^ < e Qf |w t - w^ v r d ^ = E iw* - w t st \ p dt ) ■ 

We will now show for p = 1 and p = 2 that these relations are also true in the general 
case with b ^ 0, see Subsection 2A In addition, we will get explicit expressions for 
the expected (conditional) path deviations in the case p — 2. The reason for not 
considering a general p G N is that we just want to demonstrate the phenomenon 
that the bridge versions have different sample path behavior. We also note that the 
calculations for a general pGi would be more complicated. 



2.4. Expected absolute, quadratic and conditional quadratic path devia- 
tions. First we study the L^norm E(| Wt — W^ r \) of the path deviations W t — W^ r . 

Lemma 2.5. Let (W^ r ) t £[o,T] be a Wiener bridge from to b over the time-interval 
[0,7], where beR. Then for all t G (0,T) ; 

E (\W t - Wr\) = E (\W t - W?\) > E (\W t - W?\) . 

Proof. For a normally distributed random variable Y„ c 2 with mean \i and with vari- 
ance a 2 > we have 



14 
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X 



exp 



1 (x — /i) 2 



dx 



x — jJ, 1 



a 



2,71 



exp 



1 (x — n)' 



<7 



+ /i 

By change of variables y = 



V2 



exp 



\(^) 2 and z 



1XO c 



x—fj, 



l(x- flf 



(7" 



d.r 



dx. 



we get 



E(\Y^\)=2a / — = exp(-y) dy + fi 
(2.12) = 2.-1= exp (_£) +A1 (l- 2 *(- 



/27T 



exp — z 



dz 



2a $' 



(7 



A*(2*( 



Differentiation with respect to a > 0, using = — x$>'(x) yields 

sXM = *0-^0-^(;)=*(J)>o. 

Hence E(| V^ )Cr 2 1) is a strictly increasing function in a > from which, by Subsection 



2.2 together with (2.3), we get for all < t < T, 



E (\W t - W? v |) = E - 



2t 



T 



T 



2$ 



concluding the proof. 



□ 



Next we compare expected absolute path deviations E ^ J Q T \ Wt — W t br | dt^j. Us- 
ing that integration over the time-interval [0, T] and taking expectation can be 
interchanged (as it is explained in the introduction), by Lemma 2.5, we also get 

fT 



E 



\W t - Wr\dt 



E 



r 



w t - Wr\dt > E 



\W t - W?\dt 



Using (2.12) and Proposition 2.4, it might also be possible to calculate and to 
compare expected conditional absolute path deviations given Wt = d. This task is 
more complicated, since now the mean is different for different versions of the bridge, 
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see Proposition 2.4 Instead we will now consider expected (conditional) quadratic 
path deviations which have much nicer forms. 

Next we calculate the second moments E ((Wt — W^ T ) 2 ) of the path deviations 
W t — Wt\ and also expected quadratic path deviations E ( /J \w t — W^ r ) 2 dtj. 



Theorem 2.6. Let (W^)te[o,Tl be a Wiener bridge from to b over the time- 
interval [0, T], where b £ R. Then for all i£[0, T), 



(2.13) 
(2.14) 



E ((W t - WD 2 ) = E ({W t - Wt) 2 ) = ^ + b 2 p 2 



E({W t -Wn 2 ) = a 2 (t) + b : 



rp2 ' 



where c 2 (t) is defined in Proposition 2.3. Moreover, the expected quadratic path 
deviations take the following forms: 

E Qf (W t - WD 2 dtj = E Qf (W t - W^f dtj = |(T + b 2 ), 



E 



T 



(W t - wiy dt 



T (T 



+ b< 



Proof. For a normally distributed random variable Y^ i a 2 with mean n and with vari- 
ance a 2 > we clearly have 



(2.15) 



E(Y 



Var(F MiCT2 ) + (E^)) 2 = a 2 + /i 2 . 



Hence, by Proposition 2.3, we get (2.13) and (2.14). Then, we have 



E 



(W t - WD dt 



E 



T 



(w t - wt) 2 dt 



T f 2 



T 



-(T + b 2 )dt=-(T + b 2 ) 



and, by a change of variables s = (T — t)/T and partial integration, we get 



E 



(w t -wi T ydt 



r 



At) + b 2 ^) dt 



t 2 



T 



T 2 I 

+ 2(T-t) log 
l 



T-t l2 t 2 
h b — 



dt 



T 



-T 2 + 2T 2 



slogsds + -b 2 T 



16 



mAtyAs barczy and peter kern 



■ rp2 _ rp2 



1 2 1 2 

6 3 



sds + -b 2 T 



T (T 



+ b< 



concluding the proof. 



□ 



Note that, by Theorem gj| and ggp , for all t G (0,T), 

E ((Wi - WT V ) 2 ) = E ((W, - Wf) 2 ) > a 2 (t) + b 2 p 2 = E ((W t - Wff) . 

Further, in case 6 = this shows that the expected quadratic path deviation of the 
integral representation is half of those of the anticipative version and the space-time 
transform of the bridge. This is in accordance with the typical observations from 
simulation studies as in Figure [TJ 

Next we study expected conditional quadratic path deviations. 

Theorem 2.7. Let (Wf)te[o,T] be a Wiener bridge from to b over the time- 
interval [0, T], where b G R. Then for all £ G [0, T) and d G R we have 



(2.16) 
(2.17) 



E ((W t - W t av ) 2 \W T = d) = (d- bf 



J^2 ' 



E ((W t - W?) 2 1 W T = d) = 2£ 



T-t 



T 

(T - 1) 2 
T 



(T-t) 2 T-t 
2 l _ ' log 



T 



T 



log 



T-t 



T 



(2.18) 



E ({W t - Wf) 2 \W T = d) 



( fJ . \ £ ,T-t. T-t 
+ ( {d-b)- + d^^\og^r- 

t 2 (2t-T) 2 
T f ^W** 



+ 



T T 1 



T ,T)( t ) 



Moreover, the expected conditional quadratic path deviations take the following forms: 



(2.19) E / (Wt-W^ydt 



(2.20) E / (W t -W?) 2 dt 



(2.21) eQT (W t - W t st ) 2 dt 



W T = d 
W T = d 
W T = d 



1 
3 
7 

54 

-id - b) 2 T + -b 2 T - —dbT + -T 2 . 
6 V ' 6 12 6 



(rf-6) 2 T, 
(6 - d) 2 T 



54 



— dbT + — T 2 , 
54 27 ' 
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Proof. By l\2A} , fl2.7[ ) and ( |2.15[ ), for < t < T we get ( |2.16[ ). Using that integration 
over the time-interval [0, T] and taking conditional expectation can be interchanged 
(as explained in the Introduction), we get (2.16) yields (2.19). By (2.5), (2.8) and 
(2.15), we have (2.17), and hence, by a change of variables s = (T — t)/T and partial 
integration, we get 



E 



(w t - wi r y dt 



W T = d 
t 2 



2t— + (d - 6) 2 — + 2d(d - 6)-— log — 



+ 2T^£log^ + ( d 2 -T)( T -^ 



J>2 



^2 



log 



T-t 



T 



dt 



-(d-b) 2 T + / 2T 2 (1 -s)s + 2d(d-b)T(l -s)slogs 

+ 2T 2 s 2 log s + (d 2 - T)Ts 2 (log s) 2 ] ds 

hd- b) 2 T + T 2 - \t 2 - d{d ~b)T ! sds + \d{d - b)T ! s 2 ds 
3 3 Jo 3 J 

-It 2 f s 2 ds-(d 2 -T)T f If log ads 
3 Jo Jo 3 

I(d _ & )2 T + I T 2 _ l d(d _ 6)T + _ h)T _ 2 T 2 + 2 (rf2 _ T)T ^ s2 dg 
\{d - b) 2 T + l -T 2 - ^d(d - b)T +^(d 2 - T)T 



1 

27" 



-T 



7 

54' 



d 2 T - — d&T + -6 2 T, 



which yields ( |2T20| >. Finally, by Q and fl2l5| ), we have ( |2J8| ), and hence, by 

a change of variables s = It — T, we get 

rT 

W T = d 



E (/ ( W t~ W t St ) 2dt 



-T 2 + -6 2 T - 
3 3 



(2t-T) 2 2dbt(2t-T) d 2 2 
jl 1 ^2 J^' ^ 



dt 



i T 2 + I 6 2 T _ 1 

3 3 2 



s 2 , db(s + T)s d 2 2 



+ 



^2 



rp2 S 



ds 
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3 3 2 V 3 3 2 3 



l T 2 + -b 2 T 



6 



12 



d&T + -d 2 T. 
6 



which yields (2.21) 



□ 



In what follows we give a complete comparison of the quantities (2.19), (2.20) and 
(2.21). Let b = b/y/T and d = d/y/T. Using the notation 

T 



e br :=E(Y (W t -W?) 2 dt 



W T = d 



by Theorem 2.7, we have 



e w = \(b-d) 2 T 2 



L { p_d) 2 + -b 2 --bd + -\ r 

54 v ; 54 54 27 



Hence we easily calculate 

6 av > C; r 'i )' 



l l(j b _d) 2 >-~b 2 

54 v ; 54 



-bd+-)T 2 . 
12 6 



7 rj 1 
— -bd H 

54 27 



d- —b 
22 



> 



2 

IT 



15 
22 



e av > e st <=> hb - d) 2 > H 2 - ^-bd + - 
o o 12 o 



d- -b 
4 



> \/ L + ^^) 2 



and 



6st ^ ^ir 




27 v ; 


>- 

27 


5 

108 






d--6 
8 


<v 


'> 2 


7 

~~ 2 



6rf 



54 



and (6) 2 > 



The corresponding regions are graphically illustrated in Figure [3] 



224 



SAMPLE PATH DEVIATIONS OF PROCESS BRIDGES 



19 




FIGURE 3. Regions in the (b, d)-plain for which A: e av < e; r < e s t, B: 
e; r < e av < est, C: e ir < e st < e av , and D: e av < e st < e ir . 



Finally, we remark that Theorem 2.7 justifies our simulation results in the case of 



the endpoint Wt of the Wiener sample path is close to the prescribed endpoint b 

we get 



of its bridge. Indeed, in case of d = b, by Theorem |2.7 



e(J (W t -Wr) 2 dt 



E 



E 



(w t - w?ydt 



(w t - w?y dt 



W T = d 
W T = d 
W T = d 



0. 



2 1 
27 27 

12 6 



which shows that the expected conditional quadratic path deviation of the Wiener 
process from the anticipative version of its bridge is being smaller than from the 
integral representation of the bridge or from the space-time bridge version. 



3. Path deviation of the Ornstein-Uhlenbeck process from its bridges 

3.1. Ornstein-Uhlenbeck bridge versions. Let (C/*)t>o be a one-dimensional 
Ornstein-Uhlenbeck process starting in a G K, i.e., it is the unique strong solution of 
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the SDE 

d£/" t a = qU? dt + a dW t with initial condition Uq — a 

for some g ^ and o > 0, where (W t )t>o is a standard Wiener process. It is well- 
known that the Ornstein-Uhlenbeck process has the integral representation 

U? = e qt (a + a e~ qs dW s ^ , t>0, 

which is a Gauss process with mean function E(C7") = ae qt and covariance function 
Cov(U*, U t a ) = (x 2 ^ sinh(gs) for < s < t. We also have U t a = ae qt + U°, t > 0, 
where (U^) t >o is a one-dimensional Ornstein-Uhlenbeck process starting in 0. 

We consider the following versions of the Ornstein-Uhlenbeck bridge from a to b 
over the time- interval [0, T], where a, b G M: 

1. Anticipative version 

uav = a sinh(g(T-t)) + b sinh(gt) + / _ sinh(gt) f/0^j , < t < T. 
sinh(gT) sinh(gT) \ * sinh(gT) 

Up to our knowledge this anticipative version of the Ornstein-Uhlenbeck bridge first 
appears on page 378 of Donati-Martin [TT] for a = b = and in Lemma 1 of Papiez 
and Sandison [20] for special values of q and er. It is also an easy consequence of 
Theorem 2 in Delyon and Hu [10J and of Proposition 4 in Gasbarra, Sottinen and 
Valkeila |H]. 

2. Integral representation 

stoh(,(T _ t) ) sinl^ ff /•' d ,„ if < 4 < r , 



= ^ sinh(gT) sinh(gT) J sinh(g(T — s)) 

6 if t = T. 

This integral representation of the Ornstein-Uhlenbeck bridge is the unique strong 



solution of the below given SDE (3.2), see, e.g., Barczy and Kern j3| Remark 3.9]. 
3. Space-time transform 

( sinh(g(T -t)) sinh(gt) t k{T) - K (t) 

TTSt ) a • i , ^ r b . , — W it U < t < i , 

t/£ = < smh(gi J smh(gi ) k{1 ) K (x)-«(t) 

(b if t = T, 

with the strictly increasing time-change 

, . e~ qt sinh(ot) 1 - e~ 2qt 
E9tK nit) = — = . 
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This space-time transform of the Ornstein-Uhlenbeck bridge goes back to the proof of 
Lemma 1 in Papiez and Sandison (20] and is roughly speaking a time-transformation 
by k and a rescaling with the coefficient e qt of the space-time transform representation 
(W? t )te[o,T] of the Wiener bridge from a to b over the time-interval [0,T]. 

Remark 3.1. We note that the previous versions of an Ornstein-Uhlenbeck bridge 
are in accordance with the corresponding versions of a usual standard Wiener bridge 
introduced in the introduction. By this we mean that for all T > 0, t G [0,T] and 
(7 — 1, Ut T converges to W^ r in L 2 (fi,J r , P) as q — > 0. Indeed, 

, sinh(g(T-t)) , (T - 1) cosh(g(T - 1)) T-t 

hm — = hm — = — — — , 

q^o smh(gT) q^o Tcosh(gT) T 

sinh(gt) tcoshfgt) i 

hm — = hm — = — , 

?->o smh(gTj <?->o Tcosh(gT) T 

, nt K(T)-K(t) , e" 94 - e^- 2T ) , -te-9* - (t - 2T)e^" 2T ) T-t 

hme g ^ = hm ^ — = hm — — = — , 

g^o K (T) <?->o l-e-W q ^o 2Te~ 2 i T T ' 

K(t)n{T) v / , , /t(t) 2 \ (1-e" 2 ^) 2 
lim J; ; v ; = lim /c(t) + v ; = t + lim 



</^o «(T) - /c(t) g-*> y w k(T) - «(*) ) q^o 2g(e" 2 <?* - e~ 2 i T ) 

t 2 tT 



T-t T-t 



Further, 



E(C/° - W t f = E (YV ( *~ s) ~ l ) = f*( e<i{t ~ S) ~ X ) 2 ds 



Q 2qt _ i 2 

+ -(1 - e qt ) +t -»■ as g^0, 

2q q 



and 



E i/ ^-')> dff ,-/'id dB ,;' 

/o smh(g(T - s)) J T — s 

1 fsmh(q(T-t)) T-t\ 2 , 

■ r L (( " ™ ds ^ as g -> 0, 

\smh(g(i — sJJ 1 — s J 

where the convergence follows by Lebesgue's dominated convergence theorem. □ 

In all what follows we will use the notation (f/ i br ) te [o i T] if the version of the bridge 
is not specified. 



22 MATYAS BARCZY AND PETER KERN 

First we present a lemma about a time-transformation which will be useful for 
calculating Var(f/ t a - Uf) and Cov(C/f, Uf ), < s, t < T. 

Lemma 3.2. For tne time-transformation K* T {t) := ^jz^) > * £ [0>^); w^/i K W := 
, t G R, we get zs strictly increasing and n^(t) > t for all t E [0,T). 

Proof. Since the function [0, T) 3 i K is strictly increasing and R 3 t i— >■ = 
1- ^ 2l?t is strictly increasing for every q ^ 0, we get that [0, T) 3 t (->■ ^jz^ =: 
is strictly increasing. Further, easy calculations show that 

(1) Ky(0) = and lim^ T = oo. 

(2) k* t is differentiable on [0,T), namely 

_ K '(t) K (T)( K (T)- K (t)) + K (t) K (T) K '(f) «'(t)« 2 (T) 

(«(r)-«(0) 2 "(«(r)-«W) 2 ' 

with «'(*) = e" 2<? * and hence (k^)'(0) = 1. 

(3) For the second derivative we get 

K"{t)K\T){K{T) - K{t)f + 2(«(T) - K(t))(«') 2 (t)« 2 (T) 



(k(T) - ^)) 4 
« 2 (T)(«"(0(«(T)-«(0) + 2(K') 2 (t)) 



, tG[0,T). 



("CO - ^(t)) 3 

Since «"(t)(/c(T) - «(*)) + 2(K') 2 (t) = e" 4 "* + e - 2 ^ T+t ) > we have (k*)' is 
strictly increasing. 

Altogether this shows that («5-)'(t) > («r)'(0) = 1, i G [0,T) and hence /^(i) > t for 
aUfe[0,T). □ 

Proposition 3.3. Let (£^ br )te[o,T] fre an Ornstein-Uhlenbeck bridge from a to b 
over the time-interval [0,T], where a, b e R. Tnen (f/ t br ) te [ ,T] is a Gauss process 
with mean function 

E m br) = sinh(g(T - t)) sinh(gt) < t < T 

^ J a sinh(gT) + ° sinh(gT) ' U S < ' 

and itratn covariance function 

(3.D cov ( ^, «n = £ £ s £ t < r 

Hence all the bridge versions above have the same finite- dimensional distributions. 
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Proof. For < s < t < T, we have the covariance function 



rwrrav ™ _ c ( rro _ sinh(gs) sinh(gt) 
Lov(U s ,U t )- Cov \U a ^ U T , U t ^ U T 



Gov W, ^ - ^ Cov (C / s °, B?) - ^ 00,(0?, U°) 
smh(gi ) smh(gi ) 

sinh(gg) sinh(gt) 

H • 1,2/ rriN OOV^U T , U T J 

sinn (gi ) 

2 e 9 * . sinh(gt) 2 e qT . 

a — smh(gs) a smhfgs) 

q smh(gT) q 

a 2 sinh(gs) sinh(g(T — t)) 



sinh(gT) 



and 



CoAU:, Un = o» Gov ( f dW r , /' f^g-fl dW, 

VVo smh(g(T-r)) J smh(g(T-r)) 

^ 2 r sinh(g(T - s)) sinh(g(T - £)) ^ 

7o sinh 2 (g(T - r)) 

a 2 1 

= — sinh(g(T — s)) sinh(g(T — t)) / ^— dv 

q ' Jq(T-s) sinh v 

° 2 ■ v, i frr \\ • uf (r +\\ f cos Hq( t ~ s )) cosh(gT)\ 

= — smh(g(T - s)) smh(g(T - £)) , ^ - . , , 

q \smh(g(i — s)) smh{ql ) J 

= f! s inh( 9 (T - .)) sinh( 9 (T - t) ) .y.f-f.-j'' 
g smh(g(T — s)) smh(gT) 

o~ 2 sinh(gs) sinh(g(T — t)) 

q sinh(gT) 



By Lemma 3.2 for 0<s<t<Twe get 

Cov(E/f , t/f ) = a 2 e q ^ " (T) :" (s) K{T) ~" {t) Cov ( W «,.,.,,, . !T 



... ww..,, n,^;n,^ , , , , n,(t)«(T) 

fv(i ) K{1 ) \ /t(r)-«(«) (t(T)-i«(t) 



a 2 e 



2 ?t sinh(gs) e 2qt — e 2qT a 2 sinh(gs) sinh(g(T — £)) 
e~ qT sinh(q , T) 2q q sinh(gT) 

concluding the proof. □ 

In what follows we study the continuity of the sample paths of the bridge versions. 
It follows from the definitions that all bridge versions have almost sure continuous 
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sample paths on [0,T). The (left) continuity of the trajectories at t = T is also 
obvious in case of the anticipative version, but not in case of the integral representa- 
tion and the space-time transform. The strong law of large numbers for a standard 
Wiener process (see, e.g., Problem 2.9.3 in Karatzas and Shreve |18j ) yields the de- 
sired continuity for the space-time transform. The above mentioned continuity for 
the integral representation follows from Lemma 4.5 in Barczy and Kern j3]. For the 
sake of completeness and easier lucidity, in the Appendix we formulate and prove this 
lemma in the present setting (without reference to the notations in Barczy and Kern 

my 

Hence the anticipative version [7 av , the integral representation U n and the space- 
time transform U st induce the same probability measure on (C[0, T], B(C[0, T])). 
This underlines and explains the definition of an Ornstein-Uhlenbeck bridge from a 
to b over the time-interval [0, T], by which we mean any almost surely continuous 



Gauss process having mean function and covariance function given in Proposition |3.3 
We also note that the finite dimensional distributions of the Ornstein-Uhlenbeck 
bridge versions coincide with the conditional finite dimensional distributions of the 
Ornstein-Uhlenbeck process (U^) te \o } T] (starting in a) and conditioned on {U^ = 
b}, see, e.g., Delyon and Hu [T0~| Theorem 2], Gasbarra, Sottinen and Valkeila [T4"l 
Proposition 4] or Barczy and Kern [3j Proposition 3.5]. 

3.2. Different sample path behavior of Ornstein-Uhlenbeck bridge versions. 

First we present an indicator for different sample path behavior of the Ornstein- 
Uhlenbeck bridge versions. If we consider the linear SDE 

(3.2) dU? = q (- coth(g(T - t))U* + _ dt + a dW t , 0<t<T 

with initial condition Uq t = a, then the integral representation is the unique strong 
solution of this SDE (see, e.g., Delyon and Hu [10, Proposition 3] or Barczy and 
Kern [3, Remark 3.10]) and the anticipative version and the space-time transform are 
only weak solutions. Indeed, if the anticipative version and the space-time transform 
were also strong solutions, then, by the definition of strong solution, we would get 
P([/ t av = Ul T , V t e [0,T)) = 1 and F(Uf = Uf, V t G [0,T)) = 1, which are not 
true. For example, it does not hold that Var(t/" t a — Uf v ) = Var([7 t a — Uf) for all 
* G [0,T), and Var(C/ t a - Uf) = Var(f/f - U?) for all t G [0,T), by Propositions 
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3.5 and 3^ below. We present another indicator for different sample path behavior of 
the Ornstein-Uhlenbeck bridge versions by calculating the covariances Cov(U^ r , £/") 
of the coordinates of the two-dimensional Gauss process (U^ ,U^) t ^[o,T]- Note that 
these formulas are hard to compare. 

Proposition 3.4. For all < t < T and a,b G M. we have 

a 2 sinh(g(T — £)) sinh(gt) 



cov(ur, un 



q sinh(gT) 

° 2 „-,(T-t) n ;„w„^ ^f^,^J sinh( 9 T) 



Cov(Ul r , U?) = — e' q(1 - t} sinh(g(T -t))[qt + log 



q ' \ \sinh(g(T — t)) 



V * ' * ; q K J sinh(gT) 

Proof. For the anticipative version we have 



sinh(gt) tt0 
2 e 9 * sinh(gt) 2 e qT 



C v{Ur,U?) = C a v[u°-^^U°,U t 



a — sinh(gt) - ———a — sinh(gt) 
q smh{ql ) q 

— Sin i h(gt) , (e«< sinh(gT) - e« T sinh(gt)) 
q sinh(gT) 

cr 2 sinh(g(T — £)) sinh(gt) 
g sinh(gT) 



Using Lemma 3.2 we get for the space-time transform 

Cov(0* Un = Cov {^ <T l^ {t) W^ U? 



^ <T l~ T) K{t) Gov [ [ KUt) dW„, I e^dW, 



= a e 

«(T) 

_ .2-2^ «( r ) ~ «(*) 1 ~ e " 

«(T) g 
_ a 2 sinh(g(T-t)) 
g sinh(gT) 
Finally, we get for the integral representation 
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ft e ~qs 

ctV* sinh(g(T-t)) / — — -ds 

Jo sinh(g(T - s)) 

2aV^> sinh(g(T - t)) £ j _ d, 

r 2 e -2 9 (T-t) , 

-?( T -*) s inh(g(T-i)) / — . dr 

J e -2 qT r(l - r) 



^" 9(T_t) sinh (^ - 0) (tog [^F- ) - ^ v 1 _ e _ 2gT 
^!- g( T- t )^u^^ *\{^,^{ sinh(gT) 



sinh(g(T — i)) ( gt + log I — 



^sinh(g(T — t)) 

concluding the proof. □ 



Proposition |3.4| also shows that, even though the three bridge versions have the same 
law on (C[0, T], B(C[0, T])), their joint laws together with the Ornstein-Uhlenbeck 
process U a are different. 

Our aim is to analyze the sample path deviations of the Ornstein-Uhlenbeck bridge 
versions to the original Ornstein-Uhlenbeck process (U^)t^[o,T) (starting in a) by cal- 
culating and comparing expected quadratic path deviations E ^ f Q T (U^ — U^ T ) 2 dt^j . 

Simulation studies show the same qualitative behavior of typical sample path de- 
viations of the anticipative version, the integral representation and the space-time 
transform of the Ornstein-Uhlenbeck bridge as we have for the Wiener bridge, see 
the upper row of Figure |4} Note that in general the deviation from the space- 
time transform bridge version is hard to compare with the other deviations, since 
(U^)t^[t*.T) depends on the non-visible part {U^)t<^[T,oo) of the Ornstein-Uhlenbeck 
process, where t* G (0,T) defined as follows. Due to the strict monotonicity of k* t 
and lim t -|-T = oo there is a unique t* G (0,T) such that K^(t*) = T, see the 

analysis of the time-transform in Lemma 3.2 From simulation studies we also get 
the above typical behavior is again reversed in case the endpoint Uj. of the Ornstein- 
Uhlenbeck sample path is close to the prescribed endpoint b of its bridge, namely, for 
such a sample path of the Ornstein-Uhlenbeck process the deviation from its antic- 
ipative bridge version is smaller than from its integral representation of the bridge, 
see the lower row of Figure |4} 
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FIGURE 4. Two sample paths of the Ornstein-Uhlenbeck process with a = 1 
and q = — 1 (upper row, thick line), respectively q = 2 (lower row, thick line) 
and their deviations from the anticipative version (left column), the integral 
representation (middle column), and the space-time transform (right column) 
of the Ornstein-Uhlenbeck bridge from to over the time-interval [0, 1]. 



Our aim is again to give quantitative answers to this qualitative behavior observed 
from simulation studies by studying the path deviations on [0, T): 

= (ae qT 

(3.3) U?-U? = (ae* r 
U?-U? = (ae qT 

Note that all path deviations depend only on the transformed difference (a e qT — b) of 
starting and endpoint of the bridge. Hence in the sequel without loss of generality we 
can and will assume that a = 0. For simplicity we will concentrate on calculating the 
Gauss distributions of path deviations and to compare the expected quadratic path 
deviations only. 



^ sinh(gt) + sinh(gt) ^ 



-b) 
-b) 



a 



sinh(gT) sinh(gT) 
sinh(gt) 
sinh(gT) 
sinh(gt) 
sinh(gT) 



,«(*-«) _ 



sinh(g(T — t)) 
sinh(g(T — s)) 



dW s , 



U?-ae qt 



K(T) - K(t) 

"CO 



w 



K(T)-K(t) 
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3.3. Gauss distribution of path deviations and expected quadratic path 
deviations. First we prove that path deviations have Gauss distribution. 

Proposition 3.5. Let (£/ t br )i e [o,T] be an Ornstein-Uhlenbeck bridge from to b 
over the time-interval [0,T], where i) 6 I. Then for all t G [0, T), the path 
deviation C/ t ° — £/ br is normally distributed with mean 



b _ sinh(gi) 



and with variance 
(3.4) Var(C/ t ° - U™) = a 



sinh(gT) 
,e qT sinh 2 (gt) 



q sinh(gT) ' 



qt sinh(g(T - 1)) 



+2(1 -e' 



sinh(gT) 



Proof. With a = 0, by (3.3), for every < t < T the path deviation U® — £/ br is 

s'mh(qt) 
sinh(gT) 



normally distributed with mean E(Z7° — £/ br ) = —6 sm w 9 ji and with variance 



Var(f/< ^ } " sinh 2 (gT) Var(f/T) " ° q sinh(gT) ' 

and 

Var(C/° - fn = a 2 f (e**-) - sinh(g(T ~ t}) V ds 
1 * <j io V sinh( 9 (T- S ))J 

/"* [ 2g (t-.) _ 2e 9('- s ) sinh ^( T-t )) __|_ sinh 2 (g(T - g) 
io L sinh(g(T-s)) sinh 2 (g(T - s)) 



a 2 



d.s 



1 /•* e 9(T-s) _ -q(T-2t+s) 



2q y Jo e«( T - s ) - e-<?( T - s ) 

, rt 1 

2/ — - ' 1 



+ Smh(9(r - i)) ' 3i nh^(T- S )) dS 
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a 



2q 



^ -1)- 



2(i- e - 2 ^ T -*)) f e 



Q 



v 



e q(T-t) 



-t) v 2 - 1 



dv 



. 2 1 f cosh(q(T - 1)) cosh(gT)\ 

+ smh (q(T -t))- . — -r - . 

q \smh(q(T — t)) smh(gT) / 



:i _ e _2,(r-t)) e ^ T -l 
— log 



+ sinh 2 (g(T-f)) 



q - e^T-t) _ i 

1 sinh(gt) 



q sinh(g(T — £)) sinh(gT) J ' 



2!) 



which yields (3.5). Using Lemma 3.2 we get 



Var([/° - Uf ) = Var ( U? - a e<* W_ 



Var(£7 ( 



SM3 

re(T)— re(t) 



2a ^^(T)-^ Cov( 



«CT) 



t , K(t)K(T) 

«(T)-«(t) 



= a 2 — sinh(gt) + or 2 e 29 */t(t) 

g k(T) 



o3* 



«C0 



a 2 — ( sinh(gt) + e 9 * 



1 - e- 2 " T 



-2 



e -2gt _ e -2gT 

1 - e- 2 " T 



e 9 * - 1) 



which yields (3.6). 



□ 



Note that in the proof of Proposition 3J3 below we will give different representations 
of the variances Var({7 t ° — Uf T ) calculated in Proposition 3.5 



Next we compare the second moments E((?7 t — U x 



n 2 ) 



of the path deviations 
Uf — Uf r . In view of (2.15) we have to compare the variances of path deviations for 
different versions of the Ornstein-Uhlenbeck bridge, since the mean function of path 
deviation is the same for all versions. 
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Proposition 3.6. Let (U^ v )t^\o,T] be an Ornstein-Uhlenbeck bridge from to b 
over the time-interval [0,T], where 6 G 1. Then for all te(0, T), we have 

E((U t ° - Uff) < E((C/° - C/-) 2 ), »/ 9 > 0, 
E({U? - U?Y) < E((C/ t ° - L/f ) 3 ), i/ g < 0. 

Proof. We first give different representations of Var(t/ t ° — Uf) and Var(C/ t ° — i7 t st ) 
calculated in Proposition |3.5| that are more suitable for comparison. By Proposition 



(3.7) E((0?-0*) a ) < 



3.5, we have 



Var(£/" t ° 



a 
(I 



(3.8) 



2sinh(g(T-f)) 
+ sinh(gt) ( e qt - 



sinh(gt) 
sinh(gT) 



-q(T-t) 



qt + log — 



sinh(gT) 



sinh(g(T — t)) 



q t sinh(g(T - t)) 



sinh(gT) 



2 — sinh(g(T - t)) 
Q 

r 2 n^nU^i 



sinh(gt) 
sinh(gT) 



_ e ~q(T-t) 



qt + log — 



sinh(gT) 



sinh(g(T — t)) 



+ e 



qT a sinh (qt) 
q sinh(gT) ' 



and 



(3.9) 



Var(C/ t ° - XJ X 



<7 



sinh(gt) I < 



,9* 



sinh(g(T — t)) 
sinh(gT) 



„ sinh(g(T — t)),.,,. „ 
+2 (sinh(gt) + 1 - e 9 *, 



sinh(gT) 

gT a 2 sinh 2 (gt) a 2 sinh(g(T - t))(l - cosh(g£)) 
g sinh(gT) g sinh(gT) 
The advantage of this new representation is that now the variances include the term 
g T of smh (qt) £ versions of path deviations. 
For the comparison E(C/ t ° - t/f ) 2 with E(t7 t ° - t/f v ) 2 we consider the continuous 
function /i g on [0, T] defined by 

= 2 ^ sin h fa (r- f ))(i-c 0S h M ) ( 

g smh(gT) 

Clearly, /i g (t) = if and only if t G {0,T} and further for all < t < T we have 
h q (t) < if g > and h q (t) > if g < 0. In view of ( |2.15| ) we get 

< E(([/ t ° - f/ t av ) 2 ), ifg>0, 



(3.10) 



E((C/ t ° - U? ) 2 ) 



>E((?7 t -f/r) 2 ), if g< 0. 
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sinh(gt) 
sinh(gT) 

Using that 



_ e -9(T-t) 



qt + log 



sinh(gT) 



sinh(g(T - 1)) 
smh(q(T-t)) 







< 



— cosh(gt) 
sinh(gT) 



if q < 0, 
if q > 0. 



sinh(gT) 

by log(l + x) < x, \x\ < 1, we have for all < t < T, 

sinh(gT) 



<1, te(p,T), 



sinh(gt) 

sinh(gT) 

sinh(gt) 
sinh(gT) 



r +log smh(g(T-t)); 

+ e -, ( T- t) ^ log siBh(g(r-Q) _ gf 



sinh(gT) 

g( T-t) sinh(g(T-t))-(l + gt) sinh(gT) 
sinh(gT) sinh(gT) 



sinh(gt) 

< : : rl +e 



1 



2 sinh(gT) 
1 



( e 9* - e- ? * + 1 - e- 2 ^ T -') - (1 + gf) (e* - e" 2 ^)) 



' 2sinh(gT) 

For q < it is enough to show that g q (t) > for all < t < T. Now 

g q (t) = ( e - 2 ^ - e - 2 ^ T -')) + (1 - e"*) - qt{^ - e" 2 ^) 

= (e" 2 ^*) - 1) (e"* - 1) + qte^ie- 2 ^ - l) 

< (e" 2 ^ - 1) (e"* - 1) + gte^e" 2 ^ - l) 

= (e- 2 « T - 1) (e-«* -l + qt)+ qt(e« - l) (e" 2 « T - l), 

which is obviously positive for q < and < £ < T. For g > we have to show that 
g g (t) < 2 - e qt - e~ qt for all < t < T. Now 

2 _ - e-<" - ff ,(t) = 1 - e qt + - e" 2 ^' + gte^l - e" 2 ^) 

== g q (t) 

for which g g (0) = holds and we have 

~g> q (t) = -qe* + 2ge- 2 ^ T -') - ge" 2 ^ 

+ qe«(l - e~ 2qT ) + q 2 te qt (l - e~ 2qT ) 
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for q > and < t < T, which completes the proof. Hence, by (2.15), (3.4), (3.8), 
(Kty and (pUOj), we get □ 



Moreover, by (3.7), the expected quadratic path deviations satisfy the following 
inequalities: 



[ T e((u? - c/f) 2 ) dt < fmu° t - un 2 ) At if q > 0, 

[ E((C/° - UF) 2 ) dt < [ E(([/° - Uf ) 2 ) dt if g < 0. 
■ Jo Jo 



In the next theorem we get more explicit representations of the expected quadratic 
path deviations. 



Theorem 3.7. Let (U t ^ T )te[o,T] be an Ornstein- Uhlenbeck bridge from to b 
the time-interval [0, T], where b £ R. Then we have 



over 



E 



E 



(U?-UD 2 dt 



(u?-un 2 dt 



b 2 sinh(2gT) - 2qT a 2 e qT sinh(2gT) - 2qT 

Aq sinh 2 (gT) + ~Jq 2 sinh(gT) ' 

b 2 sinh(2gT) - 2qT a 2 e qT sinh(2gT) - 2qT a 2 

a O 1 / m\ O 



Aq sinh z (gT) 



Aq 2 



^ Ta 2 cosh(gT) _ o 2 T 2 a 2 T 



q sinh(gT) 



sinh(gT) 

2rr a 2 2 

— + — e~ 2qT 
2q Aq 2 Aq 2 



<7 

(f Jo 



qT 



'l-e 



-)log^^dx, 
smh(a;j 



E 



{u»-ufydt 



a 2 e qT sinh(2gT) - 2qT + 2a 2 cosh(gT) - 1 



Aq 2 sinh(gT) 
a 2 /sinh(2gT) 



q 2 sinh(gT) 



q \ 2q 



+ T 



a 2 cosh(gT) 
2q 2 sinh(gT) 



cosh(2gT) - 1). 



Proof. For the anticipative version, by Proposition 3.5, we get 



E 



jf (U° - UD 2 dtj = J E(U? - UFf dt 



T cqT a 2 sinh 2 (gt) ^ | 
q sinh(gT) 



T b 2 sinh 2 (gt) 
sinh 2 (gT) 



dt 



e qT a 2 
q sinh(gT) 



2/ - -- u 



sinlr(gt) dt + — 



sinh (qT) 



sinh 2 (gt) dt 
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e qT a 2 /sinh(2gT) ,\ + b 2 /sinh(2gT) y 



2gsinh(gT) \ 2q J 2sinh 2 (gT) V 2g 

For the integral representation, by Proposition 3^5 and the previous calculations for 
the anticipative version, we get 

E Qf (C/° - Uff dtj = jf E([/° - L/*) 2 dt 

6 2 sinh(2gT) - 2gT 0- 2 e« T sinh(2gT) - 2qT 
~ Tq sinh 2 (gT) + ~4? sinh(gT) 

- ^ fsinhwr - i))e -«-> ( 9( + log . 3 ';";f > ) d* 

q J V. smh{q(T-t))J 

2a 2 r T 

sinh(g(T — t)) sinh(gt) dt. 



q sinh(gT) 



Here 



sinh(g(T — t)) sinh(gt) dt = J sinh(gT) cosh(gt) — cosh(gT) sinh(gt) j sinh(gt) dt 

= sinh(gT) / cosh(gt) sinh(gt) dt — cosh(gT) / sinh 2 (gt) dt 



and hence 



2a 2 



sinh(oT) [ . , , n s , cosh(oT) / sinh(2(?t) 

= — p J -hp,*) dt p - I 

sinh(gT) cosh(2gt) cosh(gT) sinh(2gt) t cosh(gT) 
4g 4g 2 

_ sinh(g(T - 2t)) t cosh(gT) 
" 4? + 2 ' 

T 

sinh(g(T — t)) sinh(gt) dt 

2a 2 ( sinh(gT) Tcosh(gT) sinh(gT)\ 
gsinh(gT) \ Aq 2 Aq J 

_o 2 Ta 2 cosh(qT) 
q 2 gsinh(gT) 
We also have, by partial integration, 

J sinh(g(T - t^-^qtdt = J 1 - 6 qt dt = | _ J te' 2 ^ dt) 

qt 2 q fte- 2 ^ T ^ f e" 2 ^"') ' 



q sinh(gT) 



4 2 V 2q J 2q 
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_ qt 2 te" 2 "( T -') e - 2 ^ T -*) 
~ ~4 4 + 8g ' 

and hence 

[' ■ ^.-rfT-t^ 2aV?T 2 T l_l_ 2gT 



/ sinh(g(T-t))e-" (T ^ ) gtdt 



Moreover, by the change of variables q(T — t) = x, we get 

2ff2 r^t^rr ^-«(T-t)^ sinh(gT) 



/ sinh(g(T-t))e- g(T -* ) log- 
7o si 



sinh(g(T — t)) 



dt 



2a 2 f qT 1 - e" 2x , sinh(gT) , 
' log ax, 



q J 2 sinh(x) 

and then we get the formula for E ( Jq(U® — Uf) 2 dt\ . We note that we are unable 
to solve the integral / Q 9T (1 — e~ 2x ) log ^j^ffi dx. 



Finally, for the space-time transform, by Proposition 3J3 and the previous calcula- 
tions for the anticipative version, we get 



E (^Vt " Uff dt) = £®(U? - Uff dt 
a 2 e qT sinh(2gT) - 2qT 2o 2 



2a 2 f T 

+ — — -/ sinh(g(T-t))(l-cosh(gt))dt 
q smh(gT) J 



Aq 2 sinh(gT) 

(1 — cosh(gT)) — / sinh(g(T — t)) cosh(gt) dt 



2a 2 / 1 ' T 



gsinh(gT) \ q 

a 2 e qT sinh(2gT) - 2qT 



+ 



Aq 2 sinh(gT) 



Here 



T 



sinh(g(T — t)) cosh(qt) dt 



sinh(gT) / cosh 2 (gt) dt — cosh(gT) / sinh(gt) cosh(qt) dt 
Jo Jo 

sinh(oT) / sinh(2oT) \ coshfoT) . , . m . 

' 1 > +T) ^-^(cosh(2gT) - 1), 



2 V 2 ? / 4 <? 

and hence we get the formula for E ^ J Q T (Uf — Uf) 2 dt) . □ 

We note that the formulas E ( f^U? - U* T ) 2 dt) are even harder to compare 
with each other than the variances in Proposition |3.5 with each other. It might also 
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be possible to calculate the Gauss conditional distribution of path deviations given 
U® = d using Theorem 2 and Problem 5 in Chapter II, §13 of Shiryaev [22j, and to 
calculate corresponding formulas for conditional quadratic path deviations. But even 
if these formulas are present, the conditional quadratic path deviations will be hard 
to compare, since they will depend on the four parameters q, b, d, T and possibly also 
on a. We renounce to give these explicit and likewise very long calculations. 



4. Appendix 

The following lemma yields almost sure (left) continuity at t = T of the integral 
representation of an Ornstein-Uhlenbeck bridge. 

Lemma 4.1. Let Tg (0,oo) be fixed and let (B s ) s > be a one-dimensional standard 
Wiener process on a filtered probability space (Q,A, (At)t>o,^) , where the filtration 
("4)t>o is the usual augmentation of the natural filtration of the Wiener process B 
(see, e.g., Karatzas and Shreve [TBI Section 5.2. A]). The process (lt)t e [o,r] defined 
by 

y . f£Sg^d*. ifte[o,n 

*' \0 ift = T, 

is a centered Gauss process with almost sure continuous paths. 

Proof. By Bauer P, Lemma 48.2], (Y t ) te \ 0} j^ is a centered Gauss process. To prove 
almost sure continuity, we follow the method of the proof of Lemma 5.6.9 in Karatzas 
and Shreve [18]. For all t G [0,T), let 

M t -= / . , , ^ ^dB s . 



sinh(g(T — s)) 

Then (Mt)te[o,T) i s a continuous, square-integrable martingale with respect to the 
filtration (At)te[o,T) and with quadratic variation 

(M) t := f , 2 )- -d S = -(coth(g(T-t))-coth(gT)), t G [0,T). 

Jo smh (q(T - s)) q 

Then lim^r(M)t = oo. By a strong law of large numbers for continuous local 
martingales, we get 

P ( lim 7— *- = 

W (M) t 
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see, e.g., Lepingle [TT[ Theoreme 1] or 3°) in Exercise 1.16 in Chapter V in Revuz 
and Yor [21 J. (We note that the above mentioned citations are about continuous local 
martingales with time- interval [0, oo), but they are also valid for continuous local 
martingales with time-interval [0, T), T G (0, oo), with appropriate modifications in 
their conditions, for such a formulation, see, e.g., Barczy and Pap [5, Theorem 3.2].) 
Then we have 

Y t = sinh(g(T - t))M t = sinh(g(T - t ))(M) t -^-, t e (0, T). 

\ M lt 

Here 

hmsinh(g(T - t))(M) t = lim - (cosh(g(T - t)) - sinh(g(T - t)) coth(gT)) = -. 
Hence we conclude P(lirn^i^ = 0) = 1. □ 
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